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ABSTRACT 

We examine the processes which could lead to the observed enhancement of 
Li and possibly other light elements (Be, B) in the companions of a number of 
X-ray novae. We conclude that one of the most promising mechanisms is the 
spallation of CNO elements on the surface of the companion induced by the 
neutron flux produced in the hot accretion flow onto the compact object. Direct 
production of the observed Li and its deposition onto the dwarf companion seem 
less likely, mainly because of the possibility of its destruction in the production 
region itself and difficulties in its deposition associated with the configuration of 
the companion's magnetic field. We discuss other potential observables of the 
above scenario. 

Subject headings: accretion — black hole physics — radiation mechanisms: 
nuclear reactions — stars: neutron — X-rays 
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1. Introduction 

High Lithium abundances have recently been detected in the secondaries of at least 
five soft X-ray transients: V404 Cyg, A0620-00, Cen X-4, Nova Mus (Martin et al. 1992, 
1994, 1996), and GS 2000+25 (Filippenko, Matheson, & Barth 1995, and Harlaftis, Home, 
& Filippenko 1996). Although the measured abundances were not extremely high, ranging 
from about 10~ 10 (for A0620-00) to about 2 x 10" 9 (for Cen X-4), they are considered 
highly unusual and surprising simply because in such late-type stars as the secondaries of 
these compact objects are, Lithium should be greatly depleted by means of nuclear burning 
within the stars. This association of enhanced Lithium abundance with the presence of 
the compact companion, raises then the possibility of the production of Lithium, either in 
the accretion flow surrounding the compact object and its deposition on the surface of the 
secondary, or its production at the latter site by processes associated with the accretion 
flow. 

In addition to this development, one of these sources (Nova Muscae) exhibited during 
its outburst a prominent gamma-ray line feature at about 450 keV (Goldwurm et al. 1992), 
which was interpreted as either a red shifted electron-positron annihilation feature, or 
possibly a nuclear deexcitation line from Helium-Helium collisions in the accretion disk 
(Martin et al. 1994; Rebolo et al. 1995; Durouchoux et al. 1996; Guessoum, Ramaty, and 
Kozlovsky 1997). In the latter reference, to enhance the line emissivity, it was assumed 
that the accreted matter was Helium-rich, such as would be the case if the companion 
star was highly evolved and had lost its Hydrogen envelope (e.g. WN Wolf-Rayet). In 
that work, it was shown that the gamma-ray line flux would be sufficiently high to be 
consistent with the 1992 observation of Nova Muscae, provided the luminosity of the 
source was high. In addition, it was pointed out that this scenario could also explain the 
high Lithium enrichment observed in such sources, provided that a sufficient fraction of 
the Lithium produced was able to escape the production region in the accretion disk and 
be deposited onto the companion's surface. Furthermore, this scenario, postulating the 
occurence of nuclear reactions in the hot accretion plasma, led to additional predictions: 
a) a delayed gamma-ray line at 478 keV should be observed upon the decay of 7 Be into 
7 Li*, knowing that 7 Be has a lifetime of 56 days, and b) that an unusual 6 Li/ 7 Li should be 
measured in these cases, the ratio varying highly (between about 10 and 1) depending on 
the temperature of the disk. 

A more recent paper (Yi & Narayan 1997) has independently proposed and investigated 
the same scenario, namely the production of Lithium in a — a collisions, in an advection 
dominated accretion flow (ADAF; Narayan & Yi 1994, 1995; Abramowicz et al. 1995; Chen 
et al. 1995), rather than a Shakura-Sunyaev disk. It was argued there that some of the 
Lithium nuclei would escape and be deposited onto the surface of the secondary, enriching 
the observed abundance. This model has the advantage of keeping the ion temperature 
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much higher than the electron temperature, thereby allowing the necessary nuclear reactions 
to proceed more efficiently The model is also convenient from the technical point of view, 
as it provides simple analytic expressions for the radial dependence of the temperature and 
the density of matter in the accretion flow as a function of the viscosity parameter a and 
the macroscopic parameters of the flow (mass of the compact object and accretion rate). 

While important in its basic setting and concepts, the treatment of Yi & Narayan 
(1997) did not take into account several processes which can drastically modify the resulting 
Li production: First and foremost, the authors did not take into account the breakup of 
both Lithium and Helium by protons (Guessoum and Gould 1989) in the hot accretion 
flow. This point is of the utmost importance since both these processes will compete 
strongly with the Helium-Helium reaction which produces 7 Li. In this respect, one should 
note that the breakup of 7 Li (by protons) into 2 cc's is an exothermic reaction, having no 
threshold energy to overcome, while at the same time, the breakup of Helium has a higher 
cross section than the a — a reaction which produces the Lithium although its threshold 
is slightly ( £ 6%) higher (see table 1); hence the same high nucleon temperatures which 
allow the production reaction to take place are more likely to destroy both the products 
and its reactant constituents. Moreover, the capture of a proton by 7 Li (giving 7 Be and a 
neutron), has a very low threshold and a very high rate. This latter reaction is normally 
in equilibrium with its inverse reaction (the capture of a neutron by a beryllium nucleus 
giving 7 Li and Hydrogen), but in the present situation, the protons in the accreting matter 
are presumably much more abundant than the neutrons, and so 7 Li is much more readily 
destroyed. The same is true also for 6 Li which is rapidly broken up by the protons. 
Secondly, the authors did not consider the use of a Helium enriched accreting matter which 
could have alleviated some of the problems outlined above: indeed, our calculations show 
that had one used a 50% H - 50% He (by mass) composition instead of their 75% H - 25% 
He, it would have resulted in an enhancement in Lithium production by a factor of 4 or 
more depending on the parameters. Finally, the analytic expressions used for the cross 
sections of the alpha-alpha production of Lithium were too approximate. The above points 
will be discussed in more quantitative manner in section 3. 

The deficiencies in the Yi & Narayan (1997) treatment mentioned above, are not 
the only impediments to materializing the proposed scenario for the enhancement of the 
Lithium abundance on the surface of the companion. There is in addition the issue of 
escape and deposition of the Li nuclei onto the companion's surface. It is difficult indeed 
to see how these nuclei, charged and relatively heavy, will overcome both the gravitational 
barrier and the magnetic fields surrounding the compact object. In thermal equilibrium, 
heavier nuclei (atoms) tend to be preferentially favored for capture rather than escape. 
One could invoke the outflows, known as jets, to explain the escape of the nuclei from the 
accretion flow (it should be noted that the mechanism of jet production is not at present 
well understood); in such a case one would have to argue, additionally, that the jet beam 
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will be intercepted by the companion star, a rather unlikely proposition given the general 
belief that the jet direction is perpendicular to the plane of the binary orbit. Last but not 
least, should the outflowing nuclei be intercepted, one way or another, by the companion 
star, it is by no means certain that they will be deposited on its surface; this depends 
strongly on the geometry of the magnetic field lines of the companion star. It is well known 
from solar flare observations that deposition of particles onto the solar surface proceeds only 
through loss cone deposition, which allows only a small fraction of the available particles 
to hit the solar surface, while galactic cosmic rays, for example, never make it to the solar 
surface, otherwise they would produce an observable flux of 7r°-decay gamma rays. 

For all the above reasons, and because we consider the correlation of enhanced Li 
abundance sites with the stellar companions of X-ray novae as significant, indicative of the 
presence of a hot (7* ~ 30 MeV) accretion flow in which nuclear processes are important, we 
examine closer a variance of the scenario proposed by Yi & Narayan (1997). Our scenario 
retains the general features of ADAF; however, rather than insisting on producing the 
required Li in the accretion flow, we consider instead the production of neutrons in the same 
setting, by the dissociation of Helium or heavier nuclei. The production of Lithium, then, 
takes place on the surface of the companion, by spallation of CNO nuclei there, induced by 
the neutrons flying from the hot accretion plasma onto the compact object. 

Neutrons are relatively easy to produce in ADAFs (either in p — a or a — a reactions), 
and a substantial fraction of them are usually produced with sufficient energy to allow their 
escape from the gravitational potential of the production site. Furthermore, being neutral, 
their escape, and also their deposition on the surface of the companion, is not affected by 
the presence of magnetic fields, and it is expected to be roughly isotropic. These properties 
allow, then, one to calculate rather reliably their escaping fraction and flux onto the surface 
of the companion, where the Lithium can be produced. Therefore, the assumption that the 
observed Li is synthesized on the surface of the companion, by neutrons produced in the 
hot accretion flow onto the compact object, greatly reduces the uncertainty of the problem, 
to effectively the viscosity parameter a of the ADAF (for a given accretion rate and mass 
of the compact object). It is worth noting that our basic assumption can be potentially 
confirmed through direct observation: If the observed Li enrichment is indeed due to 
CNONe spallation by neutrons emitted by the compact object, it should be accompanied by 
a similar enhancement in the Boron and Beryllium abundances which are also by-products 
of the same reactions. We shall come back to this point later. 

Though difficult to ascertain their presence, free neutrons are generally present in 
astrophysical settings in which processes involving nuclei of energies ^ 10 MeV are at work. 
However, their presence is of little observational consequence if they do not undergo any 
substantial interaction before they decay and/or if their energy density is small compared 
to that of the associated plasmas. On the other hand, if they do interact before decay, 
and/or their energy densities are not negligible, then their effects can be significant. Thus, 



- 6- 



the neutron 2.223 MeV proton-capture line has been observed during solar flares, while 
neutrons have also been observed directly by terrestrial neutron detectors, following such 
events (for a review see Mandzhavidze & Ramaty 1993). Furthermore, because neutron 
propagation is unaffected by the presence of magnetic fields and because, by virtue of their 
production mechanism, they can have energies and energy densities comparable to those 
of the producing plasmas, they are excellent agents of energy, momentum and angular 
momentum transport in astrophysical plasmas. Their transport properties and effects have 
been considered a number of times in the past literature. Thus, Kazanas & Ellison (1986b) 
appealed to the production of relativistic neutrons from the photodissociation of very 
energetic (10 15 eV) He nuclei, as means of transporting accretion energy from the bottom 
of an accretion flow near a neutron star to its eclipsing companion, in order to produce 
7— rays in phase with the eclipses' ingress and egress. Giovanoni & Kazanas (1990) and 
Contopoulos & Kazanas (1995) considered relativistic neutrons as a means of powering the 
relativistic jets observed in Active Galactic Nuclei (AGN). Atoyan (1992) has considered the 
effects of neutrons produced near the compact object in AGN on the accretion flows which 
power them, while Guessoum & Kazanas (1990) have explored their potential importance 
as the viscous agents of accretion flows onto compact objects. 

The present treatment is much more limited in its scope: It assumes the presence of an 
ADAF similar to that of Yi & Narayan (1997) in order to model the flow onto the compact 
member of the X-ray binary (nova), and examines in detail its consequences within the 
scenario involving the production of neutrons outlined above. In §2 we briefly describe 
the background flow, its scalings and range of the parameters involved. We also discuss 
and give a list of the most prominent reactions associated with the interaction of nuclei 
in the flow and in particular the production of neutrons. In §3 we discuss the details of 
the neutron production, as well as of several other nucleids, in the ADAF along with the 
expected gamma-ray line fluxes. In §4 we discuss the production of Li on the surface of the 
companion given our computed flux of neutrons and in §5 the results are summarized and 
discussed. 



2. The Accretion System 

2.1. Accretion Model and Parameters 

The nuclear reactions relevant to our problem and to the production of gamma rays 
necessitate energies of about 20 MeV per nucleus; therefore, these processes can potentially 
take place only in the inner regions of an accretion flow, provided additionally, that the ion 
temperatures could reach values of at least several MeV. The latter condition, therefore, 
requires that the emitting plasma be optically thin, since the standard optically thick - 
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geometrically thin Shakura-Sunyaev (1973) disk is much too cool for the required nuclear 
reactions to take place. The Advection Dominated Accretion Flows (ADAF) of Narayan & 
Yi (1994), which allow for ions with roughly virial temperatures, provide a natural setting 
for our calculations. We therefore consider the regions around the compact object between 
about 100 Rs (where Rs is the Schwarzschild radius of the object, Rs = 2GM/c 2 ) and 
either the surface of the neutron star or the horizon of the black hole, depending on the 
nature of the compact object. In order for our treatment to remain generally applicable 
to both possibilities, we assume a mass of 1 M Q for the compact object, unless a specific 
value is considered. Because in the present treatment we are only interested in the nuclear 
reaction rates and not in the entire emitted spectrum, the value of the electron temperature 
does not affect our results, which thus depend solely on the value of the ion temperature TV 

In the present investigation, and for the reasons discussed above, we use the ADAF 
prescription to determine the dynamics and thermal properties of the accreting fluid; hence 
the ion temperature Tj is given by the expression (Yi & Narayan 1997) 

Ti = 1.11 x 10 12 r _1 K = 95.68 r" 1 MeV ~ 0.1 r' 1 m p c 2 (1) 

where r is the dimensionless radial variable r = R/Rs- 

Within ADAF, the density in the hot plasma is given by an analytic expression derived 
from mass flux conservation; this has a form very similar to that for spherical accretion and 
it is given by (Yi & Narayan 1997) 

n = 6.4 x 10 18 a" 1 m' 1 m r" 3/2 cm" 3 (2) 

where m is the dimensionless mass measured in units of a solar mass (m = M/M ), rh 
is the dimensionless accretion rate measured in units of the Eddington accretion rate 
m = M/M E dd = M c 2 /L E dd = M/1.4 x 10 17 m g/s, and where, following Yi and Narayan, 
we have taken (3 (the ratio of gas pressure to total pressure) to equal 1/2. 

The largest, by far, part of the kinetic energy associated with the inflow is carried 
by the hot protons which provide most of the accretion inertia. Their energy is given off 
to the electrons via Coulomb collisions (or some other type of plasma interaction) and 
then transferred efficiently to photons to produce the observed radiation. In this scenario, 
therefore, the observed luminosity depends on the rate of energy transfer from the hot protons 
onto the electrons (provided that the associated time scale is longer than the local free-fall 
time). If a v is the rate of the relevant reactions, the total luminosity of the system is given 
by L ~ / n(r) 2 a v kT { dV where dV = 4irR 2 dR = An R 3 S r 2 dr ~ 3.4 x 10 17 m 3 r 2 dr cm 3 is 
the volume element in the accretion flow. The resulting luminosity, then, will be 



L ~ 4.4 x 10 38 kT \ aV a~ 2 m m 2 
rripC 2 <7t c 
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~ 4.4 x 10 37 a 2 mm 2 erg sec 1 (3) 

where in the last step the value of k% was substituted using equation (1), with r in denoting 
the innermost radius for which the scaling of equation (1) holds and ctt, the Thomson 
cross section, was introduced to provide a normalization of the reaction rate a v to ot c. 
Assuming m to be in the range of 1 - 10, the entire uncertainty of the problem is in fact 
reduced to the value of a (following Narayan & Yi 1994 we consider a = 0.3, 0.1) and to a 
choice for the accretion rate m or equivalently M, which can be obtained from the observed 
total luminosity of the system; in the present note we will consider that to be in the range 
of 10~ 12 M Q /yr and 10~ 8 M Q /yr. 

Finally we note that the system is taken to be in steady-state, whereby the accretion 
and in-fall are constant in time, and thus the neutron production and gamma-ray emission 
are steady. 



2.2. The Nuclear Reactions 

Assuming an initial composition for the accreting plasma, usually some proportion 
of Hydrogen and Helium, we set up a network of nuclear reactions which the ions (the 
primary protons and alphas as well as the secondary species) can undergo, including those 
leading to de-excitation gamma-ray photon production; we then let the plasma evolve over 
the dynamical time scale tj with its temperature following values given by Equation (1) as 
it sinks in the gravitational well. The main reactions in our simulations are displayed in 
Table 1 along with the center-of-mass threshold energies necessary for these reactions to 
proceed. These energies serve as a rough indicator of the relative importance of the various 
channels (since the cross sections are approximately of the same magnitude). 

The variations in the abundances of the nuclear species are treated through a reaction 
rate formulation, using an expression for binary processes that is correct for non-relativistic 
energies: 

where i and j are the interacting nuclei in a particular reaction, u{vij) its total cross section, 
v y is the relative velocity of the two species, 8^ is the Kronecker delta and drii , drij are the 
distribution functions of the particles in velocity space. We assume that the timescale for 
thermalization among the ions is much shorter than that for the nuclear reactions or the 
Coulomb collisions, so that the particles' distribution functions are always Maxwellian (see 
Guessoum and Gould 1989 for a discussion of this point). We use experimentally- measured 
cross sections to obtain the reaction rate of each nuclear reaction as a function of the ion 
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temperature (references for the cross sections are given in Guessoum and Gould 1989; 
Guessoum and Kazanas 1990; see also Glagola et al. 1982). We then use a simple numerical 
scheme of explicit finite-differencing to follow the abundances of the various species as well 
as the photons produced, as the ions sink in the gravitational well of the compact object 
with the temperature changing according to equation (1), until they are removed from the 
plasma at t — t d . This calculation is performed for various values of the model parameters: 
M = 1(T 12 , 10" 11 , 1(T 10 , 3 x 1(T 10 , 1(T 9 , 3 x 1(T 9 and 1(T 8 M /yr; a = 0.3 and 0.1; 
M = 1, and 10; a Helium fraction Y = 10%, 30%, 90%, etc. 

Concerning the gamma-ray lines, since we assume the plasma not to contain any nuclei 
heavier than Helium, then the only lines of interest in our case are the following: 

• 0.429 MeV resulting from the de-excitation of 7 Be*, which is produced mostly in a — a 
collisions; 

• 0.478 MeV resulting from the de-excitation of 7 Li*, which is produced mostly in a — a 
collisions, plus a contribution from 7 Be, which, when it decays, produces 7 Li* 10% of 
the time; 

• 3.561 MeV resulting from the de-excitation of 6 Li*, which is produced mostly in a — a 
collisions; 

• 2.224 MeV resulting from the capture of neutrons on protons, where the neutrons and 
some of the protons are produced in the breakup of alphas. 



3. Production of Neutrons, Nuclear Species, and Gamma-Rays 

3.1. Neutrons 

Once the reaction rates for the network of nuclear reactions have all been calculated, 
the numerical scheme is run for each set of parameters, assuming that the temperature 
profile is constant in time and given by equation (1). The normalization of the neutron 
and nuclear line photon production rates are then given by integrating the reaction rates 
(computed as discussed above) multiplied by the square of the local plasma density, over 
the volume of the flow within which the plasma temperatures are above the reaction energy 
thresholds, indicated in the previous section. This prescription leads to the following value 
for the total rate of production of a given species (be that an ion or a photon) 

Rt ot ~3x 10 41 a~ 2 m m 2 lnr M sec -1 (5) 
a T c 
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where r M is the maximum radius within which the relevant reaction is energetically allowed. 

One should note that the radial density dependence of ADAF guarantees a rather 
weak explicit dependence of the rates on the outer edge of the reaction volume; most of 
the dependence is implicit through the reaction rate a v. It also indicates that the total 
reaction rate is proportional to the mass of the compact object and to the square of the 
normalized accretion rate m (rather than being proportional to it; this is a property of 
ADAF and generally any accretion flow in which the emitted luminosity is the result of 
binary interactions within the flow; see e.g. Kazanas & Ellison (1986a)), or in view of 
equation (3), proportional to the bolometric luminosity L. 

Tables 2a and 2b present the results of these calculations, showing the rates of neutron 
production and escape for various values of the viscosity parameter a and the mass accretion 
rate M. All the results shown correspond to M — 1M©, and R pz 100 R,s ~ 3 x 10 7 cm. 
Table 2a corresponds to a solar initial composition (i.e. 90% H and 10% He by number) of 
the plasma, while Table 2b corresponds to a 90% He, 10% H composition, appropriate for a 
He-rich companion. 

It is interesting to note that the resulting neutron fluxes depend quite strongly on the 
3 parameters considered, namely a, M, and the initial composition of the plasma. For the 
smallest values of M, the neutron flux scales roughly as its square (and also as a~ 2 ), as 
demanded by the corresponding scaling of the emission measure. However, for higher values 
of this parameter the flux saturates as all the available Helium nuclei break up and release 
their neutrons. These saturation values are in fact quite high when one considers that they 
correspond to the fluxes of the escaping neutrons. 



3.2. Li, Be, 3 He 

The numerical scheme described above can also keep track of the variation of the 
abundances of all the nuclear species of relevance, such as 7 Li, 7 Be, 6 Li, 3 He, etc. One can 
compute the final abundances of these isotopes for each set of parameters and compare them 
to each other or to neutron abundances in order to assess the efficiency of the production of 
each nucleus compared to that of neutrons. 

Our results are presented in Tables 3a and 3b where we show the relative abundances 
of neutrons and 7 Li, for various values of the viscosity parameter a and the mass accretion 
rate M. As for Tables 2a and 2b the results shown here all correspond to M — 1M Q and 
R fa 100 Rs ~ 3 x 10 7 cm. Again Table 3a is for a solar initial composition (i.e. 90% H and 
10% He by number) of the plasma; Table 3b is for a 90% He and 10% H composition. 

It is interesting to note that final neutron abundances are much higher than those of 
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Lithium, although the values given in these tables do not take into account the eventual 
escape of the neutrons, and therefore these numbers are useful only for comparative 
purposes. 

We have also estimated the error made by using the approximate cross sections 
of Yi & Narayan (1997) for the production of Li by He- He collisions, instead of the 
experimentally measured cross sections we have used (Glagola et al. 1982). Their expression 
for the cross section overestimates the actual data by about a factor of 2 at low energies 
(8.5MeV < E < AOMeV) and by a factor of about 10 at higher energies (E > 50MeV). 
This then translates into an overestimate of the Li abundance achieved in the disk by as 
much as an order of magnitude. 

Moreover, we have calculated the effect of Lithium and Helium breakup in the flow on 
the final Li abundance by running our numerical code cases with or without the destruction 
reactions. We find that for low accretion rates (m £j 10 -3 ) the difference is small, but for 
substantial rates (m £ 10~ 3 ) neglecting the breakup results in an overestimation of the 
abundance of 7 Li by a factor of up to 30 at the highest end of the accretion rate range. 
This sensitivity on the accretion rate is to be expected: Li production, and generally 
all binary reactions in an ADAF, proceed vigorously when the inverse reaction rate 
becomes comparable or shorter to that of free fall, i.e. for large m. However, while the Li 
production rate depends on the square of the He density (and thus is proportional to m 2 ), 
its destruction rate is proportional to the product of the Li and H densities, i.e. to m 3 , 
leading to the computed much reduced final Li abundance at high m. 



3.3. Gamma-Ray Fluxes 

As mentioned in the introduction, one can also calculate, through the same numerical 
program, the fluxes of the nuclear gamma-ray lines emitted in the accretion disk. For 
this, one needs only the cross section for the excitation of 7 Li, 7 Be, and 6 Li either in 
a — a collisions or in collisions of 7 Li, 7 Be, and 6 Li with protons; for lack of experimentally 
measured cross sections, we disregard the latter possibility. 

Table 4 shows the sum of fluxes of the 0.478 keV and 0.429 keV lines, since they 
combine to form a single broad feature (their Doppler widths being of the order of ~ 
100 - 300 keV). The fluxes are again shown for various values of the viscosity parameter 
a and the mass accretion rate M. However, because the fluxes are extremely low for 
normal abundances, we only show here the results corresponding to a 90% He and 10% H 
composition. The other parameters retain the same values. 



We note that the fluxes obtained here are much lower than those in Guessoum, 
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Ramaty, & Kozlovsky (1997), first because the accretion rates assumed here are very low 
(10~ 12 , 1CT 11 , 1CT 10 , etc.), and second because the density of matter in the ADAF model 
used in the present investigation is much lower than that used in the above reference. 



4. Production of Li and Other Elements on the Companion Star 

The neutrons having been produced in the accretion disk, they can escape rather 
easily. The escape is handled the same way it was done in Aharonian & Sunyaev (1984) and 
Guessoum & Kazanas (1990), that is, one simply needs to evaluate that fraction by counting 
the neutrons that have enough energy to overcome the gravitational potential energy, while 
taking into account the possibility of elastic neutron-nucleon collisions, which will impede 
the neutron escape. For simplicity we assume the neutrons to have a thermal distribution 
of temperature equal to that of the protons and the a's, the temperature being that of the 
shell were they are produced (i.e. between 10 and 100 MeV); we believe this assumption to 
be reasonable, especially at high accretion rates and/or temperatures less than say 50 MeV, 
when the time scale for thermalization becomes small compared to the free fall time scale. 
Assuming further that the neutrons are produced and fly out isotropically, one also needs 
to average over angles. Escape fractions are found to range between 10 - 30 %. 

The neutrons could in principle also decay in flight, but the sources we are interested 
in have companion stars sufficiently close to their accreting compact objects (the orbital 
periods are about a few hours, indicating separations of order of 10 11 cm) that the 
beta-decay of neutrons is relatively unimportant. 

The above considerations indicate that one can calculate in a straight forward way the 
flux of neutrons impinging on the surface of the star companion to the compact object. 
The neutrons will then undergo two types of processes: inelastic collisions which degrade 
their energies leading to their eventual capture by hydrogen, and spallation of the heavier 
elements (C, N, O, Ne, etc.). To compute the amount of Li, Be, B produced by the 
latter reactions one needs to know the relative abundances of the heavier elements in the 
atmosphere of the companion star, which can in fact vary significantly depending on its 
evolution stage, its history, etc. and though we have raised the possibility of a He-rich star 
atmosphere, in the following calculations we will limit ourselves to the solar composition 
case. 

The nuclear spallation reactions which the neutrons can induce are presented in Table 
5, along with the threshold energies for each channel. As previously noted, these spallation 
reactions will not only produce Lithium, but also Beryllium and Boron. It is also possible 
in principle to break the latter species into Lithium and Helium, but this will take a much 
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longer time and will be presently ignored. 

Finally, in order to compute the resulting light element abundances one needs the cross 
sections for these specific reactions and channels. Unfortunately, experimentally-measured 
values for these are not available in the literature, so one is forced to use estimates. In order 
to remain on the conservative side, we assume the neutron-induced spallation reactions 
to have half the proton-induced spallation reactions cross sections; the latter have been 
compiled by Ramaty, Kozlovsky, Lingenfelter, and Reeves (1997). These cross sections 
typically rise very rapidly with energy near threshold, then take constant values usually 
between 10 and 30 mb. 

We can now compute the fraction of CNONe nuclei which will be converted into 
Li when exposed to a given neutron flux. Assuming a value R n for the total neutron 
production rate, D for the separation between the compact object and the companion, 
Xc,N,o,Ne for the abundance of the corresponding C, N, O or Ne on its surface, and a sp the 
spallation reaction cross section, the fraction / of these elements which will be spalled after 
exposure to a neutron flux corresponding to these values of the relevant parameters for a 
time t will be 

/ = 7~p2 a sp T Xc,N,0,Ne • (6) 

Assuming as fiducial values for the neutron production rate the value R 39 = (i? n /10 39 ) 
neutrons/second, for the time scale r = (r/1 year), and for the distance (Dn = D/10 11 
cm) we obtain Q 



/ ~ 10- 5 -% X_3 • (7) 

The above estimate indicates that in order to enhance Li to the observed values one 
needs to spall roughly 10~ 7 — 10~ 6 of the available C, N, O, Ne nuclei (the value of X in 
the above equation has been normalized to the sum of the C, N, O and Ne abundances). It 
appears from the above figures that a short exposure to the above flux over a short flare 
lasting roughly a month would suffice to produce the observed enhancement. Of course, 
the neutron fluxes are expected to be as high as the fiducial value only during times of 
outburst, which last for several months and occur roughly once a decade. 

Yi & Narayan (1997) have discussed in some detail the question of depletion of 
Lithium from the atmosphere of the companion star once it has been deposited/produced 
there. The time scale they have used for this is ~ 10 7 yrs (they quote estimates from the 



4 The semi major axis of most of the systems discussed in the introduction is a few times 
this value, while it is an order of magnitude larger for V404 Cyg 
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literature ranging between 10 6 and 10 8 yrs). Clearly this time scale is much too long to 
alter the conclusions of our model given above. Of course, the more relevant time scale 
might simply be the convection time scale of the companion star's convective envelope, 
since that is the time scale at which material is mixed. That mixing time scale is of order 
r r ~ R*/v ~ 10 6 i?ii/t>5 s, where the velocity of convection has been taken to be 1 km/sec. 
That means however that the entire convective envelope will be mixed on similar time scales 
(10 6 s), and thus irradiation over the duration of an X-ray outburst (several months) will 
result in an enhancement in the Li abundance at about the observed levels throughout the 
entire convective envelope. Depletion (destruction) of the Lithium happens only through 
convective overshoot, which allows matter from the envelope to penetrate to deeper, hotter 
regions where the Lithium can be destroyed (Sweigart, private communication), to which 
the time scale quoted by Yi & Narayan (1997) refers. So, if these outbursts last only a 
few months (and not millions of years), convective mixing will only homogenize a given Li 
overabundance over the volume of the envelope, while convective overshoot will not have 
time to get this Li destroyed. 

In our view, a more stringent constraint on the Li abundance built up comes from the 
fact that the surface layers of the companion are gradually removed to be accreted onto 
the compact object. It is hence useful to estimate the time scale it takes for the exposed 
material on the companion's surface to be removed and accreted onto the compact object. If 
n(h) is the matter density as a function of depth from the stellar surface, the neutrons will 
penetrate to a depth at which a J n(h) dh ~ 1, with a ~ 10 mb = 1CT 25 cm~ 2 . Assuming 
a radius R = 10 u Rn cm for the companion, the mass of the surface layer affected by the 
neutrons, within which one expects an enhanced light element abundance, is 



This amount of mass would require about a year to be removed, if accreted at the 
Eddington rate onto a compact object of m — 1. It thus becomes apparent that, within the 
present model, there could very well be a variability in the observed abundances, with the 
highest abundances most likely to be observed shortly after a major outburst, before the 
processed matter is accreted back onto the compact object. One might argue that variation 
in the observed abundance due to the removal of these surface layers would not take place 
because the removed Li would be regenerated by the neutrons produced in the hot ADAF. 
However, this is not so, because while the removal rate is proportional to rh the neutron 
production in an ADAF is proportional to m 2 . Thus, following the vestiges of the present 
model, one concludes that monitoring the potential variability of the Li, Be, B abundances 
in X-ray novae systems could provide a long term probe of the accretion mode onto the 
compact objects. 




(8) 
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5. Summary and Conclusions 

In the present note we have reexamined in greater detail the scenario put forward by Yi 
& Narayan (1997) for the production of the Li overabundance, observed in the companions 
of the X-ray novae, in the associated hot accretion flow around the compact object. While 
we think that their general notion of a hot ADAF as the fundamental cause of the observed 
Li enhancement is basically correct and in agreement with the apparent association of 
these objects with the X-ray novae companions, we have opted for a variation of the above 
scenario. In particular, we have proposed that the enhancement of the Li abundance takes 
place on the surface of the late type, low mass companions of these systems by the high flux 
of neutrons produced in the hot ADAF. We believe this to be a more promising scenario as 
it does not suffer from the problems associated with the destruction of the Li in the ADAF 
itself, or the uncertainties of its transport from the accretion flow onto the surface of the 
companion star. 

We have also suggested and investigated the possibility that Li abundance enrichment 
would, in such scenarios, be enhanced if the accreted matter were Helium-rich, such as 
would be the case if the secondary were a highly evolved star (like a Wolf-Rayet) with an 
atmosphere consistently in large part of Helium and/or Carbon and Oxygen. This idea 
is indeed a speculative one, and there is no direct observational evidence that any one of 
the secondaries of the soft X-ray transients is or resembles a WR star; however, optical 
studies of Cen X-4 have led Chevalier et al. (1989) to suggest that the secondary in this 
case is probably a stripped giant, i.e. one with a rich-He atmosphere. It is interesting to 
note that of the 5 SXT's observed to have high-Li abundances, Cen X-4 has the highest 
such enrichment. We should point out here that our calculations of the production of 
Li at the companion star (section 4) assumes a normal, solar composition of the star 
(X c> N,o,Ne ~ 1CT 3 ). This possibility of Helium enrichment in the accreted matter could 
reveal itself in future observations. 

Perhaps one of the most interesting consequences of our scenario is the prediction of the 
presence of other light elements, namely Be and B, in the same sites where the Li is found 
(see Table 5). Unfortunately, the most prominent lines of these elements are in the UV and 
therefore not accessible with ground telescopes. One of us, (DK) along with Chris Shrader, 
have obtained spectra of these objects in the UV using GHRS on HST, in search of the Be 
and B lines, considering the Li line to be emitted from the accretion disk itself. However, 
at the time of the observation the UV disk emission from these objects was negligible and 
in the meantime it became apparent that the Li line emission is actually associated with 
the companion to the compact object. Since the companions to these objects are usually 
low luminosity, late type stars, the UV flux in our observations was too small to provide 
the high S/N measurements required for such detections. A more favorable prospect for 
the detection of the Be and B lines is to search for them in High Mass X-ray Binaries, like 
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Cyg X-l which are bright in the UV and which are thought to contain ADAFs favoring the 
production of neutrons, similar to those discussed above (Yi & Narayan 1997). 

An additional conclusion of the present model which has potential observational 
consequences is the possibility of time variability of the observed light element (Li, Be, B) 
abundances. One expects these abundances to be highest shortly after a major outburst, 
thus making these epochs the most favorable for their detection. We hope that this paper 
will stimulate further searches in the UV for the transitions associated with Be and B, also 
for less prominent transitions in the optical domain, but with the larger telescopes presently 
available, as well as the potential time variability of the Li abundance. 

The prominent neutron fluxes indicated by our model raise, in addition, the issue of 
production of the 2.223 MeV gamma-ray line from the capture of neutrons onto protons 
and the possibility of its detection. This line has been observed in solar flares, the only 
astrophysical site where neutron production has been observationally verified. Assuming 
that each neutron results in a 2.223 MeV photon the flux computed in the previous sections 
(Tables 2a and 2b) can be used to estimate the associated 2.223 MeV photon production 
rate by simply multiplying the associated production rate of neutrons by the solid angle 
subtended by the companion of the compact object. It is estimated that the latter is of order 
~ 1CT 3 , leading to a flux of only 10 36 — 10 37 photons per second for an compact object of ~ 1 
solar mass; for a source at a distance of 1 kpc this corresponds to a flux of 10 -8 — 10 -7 cm" 2 
sec -1 , too low to be observable with near future instrumentation (INTEGRAL). However, 
for a compact object of larger size, the flux of the 2.223 MeV line could be observable; for 
instance, assuming m — 10, a — 0.3 and m — 1, equations (3) and (5) yield a luminosity 
L ~ io 38 - 5 erg sec -1 and a neutron production rate of R ~ 10 42 - 5 sec -1 . Assuming the 
same solid angle of interception by the companion as above (~ 10~ 3 ), this figure yields the 
respectable flux of ~ 3 x 10~ 5 photons cm -2 sec -1 , suggesting that an observation of this 
line may not be totally improbable with the instruments aboard the INTEGRAL mission. 

A potential enhancement of the 2.223 MeV flux beyond the figures given above would 
be possible, if at larger distances the density profile of the accretion flow were to become 
less steep than the r -3 / 2 dependence of ADAF and result to a column density sufficiently 
high to allow a fraction of the escaping neutrons to be deposited, cooled and captured 
within the flow in collisions with the accreting gas. Under such circumstances, most of the 
escaping neutrons could be used in the production of 2.223 MeV photons, increasing the 
probability of detecting this nuclear line. 

We would like to thank Reuven Ramaty, Chris Shrader, Benzion Kozlovsky, and Allen 
Sweigart for a number of helpful discussions and comments pertaining to the present work, 
as well as the referee for several useful suggestions. NG would also like to thank Reuven 
Ramaty for the hospitality at GSFC where most of this was performed. 
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TABLE 1 

Nuclear Reactions in the Disk and their Threshold Energies 



Reaction Number Nuclear Reactions Threshold Energies [MeV] 



1.1 p+ a — > 3 He + d 18.35 

1.2 — > 3 H + p + p 19.81 

1.3 — > 3 Hc + p + n 20.58 

2.1a a + a — > 7 Be + n 18.99 

2.1b — > 7 Be* + n 19.42 

2.2a — > 7 Li + p 17.35 

2.2b — > 7 Li* + p 17.83 

2.3a — > 6 Li + p + n 24.60 

2.3b — > 6 Li* + p + n 28.16 

2.4 — > a + 3 He + n 20.58 

2.5 — >a + 3 H + p 19.81 

2.6 — > a + d + d 23.85 

3.1 a + n — > 3 H + d 17.59 

3.2 — > 3 H + p + n 19.81 

3.3 — > 3 He + n + n 20.58 

4.1 p + 7 Li — > a + a -17.35 

4.2 — > 7 Be + n 1.65 

5.1 n + 7 Be — > a + a -18.99 

5.2 — > 7 Li + p -1.65 

6.1 p + 6 Li — > 4 He + 3 He -4.02 

6.2 n + 6 Li — > 4 He + 3 H -4.78 

7.1 p + 3 He — > d + p + p 5.49 

7.2 — > 3p + n 7.72 

8.1a a + 3 He — > 6 Li + p 4.02 

8.1b — > 6 Li* + p 7.58 

8.2 — > 3 He+ 3 He+n 19.81 

8.3 — > 3 He+ 3 H+p 20.58 
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8.4 

9.1 
9.2 

10.1a 
10.1b 
10.2 
10.3 
10.4 

11.1 
11.2 

12.1 

13.1 



p + 3 H 
a + 3 H 



a + d 

n + 3 He 
p + n 



3H e + d + d 


23.85 


3 He + n 


0.76 


d + d 


4.03 


6 Li + n 


4.78 


6 Li* + n 


8.35 


3 He+ 3 H+n 


20.58 


3 H+ 3 H+p 


19.81 


3 H+ d + d 


23.85 


3 He + 3 H 


14.32 


a + p + n 


2.22 


3 H + p 


-0.76 


d + 7 


-2.22 
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TABLE 2a 

FLUXES OF NEUTRONS (for an initial 90% H and 10% He Composition) 

(neutrons s _1 ) 



M (M /yr) 


a = 0.3 


a = 0.1 


io- 12 


1.5 x 10 35 


3.5 x 10 36 


io- n 


4.0 x 10 37 


3.0 x 10 38 


io- 10 


3.3 x 10 39 


1.8 x 10 39 


3 x IO" 10 


7.6 x 10 39 


4.9 x 10 39 


io- 9 


4.0 x 10 39 


1.6 x 10 40 


3 x IO" 9 


1.2 x 10 39 


7.8 x 10 39 


io- 8 


8.7 x 10 38 


3.2 x 10 39 




TABLE 2b 


FLUXES OF NEUTRONS (for an initial 10% H and 90% He Composition) 




(neutrons 




M (M /yr) 


a = 0.3 


a = 0.1 


io- 12 


6.6 x 10 36 


5.4 x 10 37 


io- 11 


5.8 x 10 38 


1.8 x 10 39 


io- 10 


1.8 x 10 40 


2.2 x 10 40 


3 x 10~ 10 


7.1 x 10 40 


6.0 x 10 40 


io- 9 


2.1 x 10 41 


2.0 x 10 41 


3 x IO" 9 


9.4 x 10 40 


9.7 x 10 40 


io- 8 


4.1 x 10 40 


5.2 x 10 40 
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TABLE 3a 

Final 7 Li / Neutron Abundances in the Disk 

(for an initial 90% H and 10% He Composition) 



M (Mq/jt) a = 0.3 a = 0.1 

10~ 12 8.6 x 10~ 6 / 4.1 x 10" 4 6.0 x 10~ 5 / 1.0 x 10~ 2 

10- 11 5.8 x 10~ 5 / 1.1 x lO- 2 9.1 x 10~ 5 / 7.9 x 10~ 2 

10~ 10 8.8 x 10~ 5 / 7.7 x 10~ 2 1.4 x 10~ 4 / 2.3 x 10" 1 

3 x 10~ 10 8.7 x 10~ 5 / 7.6 x 10~ 2 1.3 x 10~ 4 / 2.1 x 10" 1 

10~ 9 5.3 x 10~ 5 / 4.3 x 10~ 2 7.2 x 10~ 5 / 1.4 x 10" 1 

3 x 10- 9 4.0 x 10~ 5 / 3.1 x 10~ 2 5.9 x 10~ 5 / 1.1 x 10" 1 

10~ 8 2.9 x 10~ 5 / 2.9 x 10~ 2 4.2 x 10~ 5 / 8.5 x 10~ 2 



TABLE 3b 

Final 7 Li / Neutron Abundances in the Disk 

(for an initial 10% H and 90% He Composition) 



M {Mq/jt) a = 0.3 a = 0.1 

10~ 12 7.2 x 10~ 4 / 1.9 x 10~ 2 4.6 x HT 3 / 1.2 x 10" 1 

lO" 11 5.0 x 10~ 3 / 1.4 x lO" 1 1.7 x 10~ 2 / 2.7 x 10" 1 

10~ 10 1.9 x 10~ 2 / 2.8 x lO" 1 3.2 x 10~ 2 / 3.9 x 10" 1 

3 x 10~ 10 1.9 x 10- 2 / 2.8 x lO" 1 3.1 x 10~ 2 / 3.7 x 10" 1 

10- 9 1.1 x 10~ 2 / 1.7 x 10- 1 1.7 x 10- 2 / 1.9 x 10- 1 

3 x 10"9 8.3 x 10~ 3 / 1.3 x 10" 1 1.4 x 10~ 2 / 1.5 x 10" 1 

10~ 8 6.1 x 10~ 3 / 9.8 x 10~ 2 9.9 x 10~ 3 / 1.0 x 10" 1 
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TABLE 4 

FLUXES OF 0.478 + 0.429 MeV LINES 
(photons cm 2 s _1 ) 



M (M /yr) a = 0.3 a = 0.1 



io- 12 


7.6 x 10" 


ii 


5.3 x 10" 


-10 


io- n 


5.7 x 10" 


-9 


1.8 x 10- 


-8 


io- 10 


1.8 x 10- 


-7 


2.5 x 10- 


-7 


3 x 10~ 10 


6.7 x 10" 


-7 


7.7 x 10" 


-7 


io- 9 


2.5 x 10- 


-6 


2.6 x 10" 


-6 


3 x IO" 9 


5.8 x 10- 


-6 


6.2 x 10- 


-6 


io- 8 


3.7 x lO- 


-6 


4.9 x 10- 


-6 
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TABLE 5 

Neutron-Induced Spallation Reactions of C, N, O, and Ne 

Nuclear Reactions Threshold Energies [MeV] 



n + 12 C — > n + 3a 7.28 

— > a + 9 Be 5.70 

— > 3 Hc + 10 Be 19.47 

— > 3 H + 10 B 18.93 

— > 7 Li +a + 2 H 22.40 

n + 14 N — > a+ n B 0.16 

— > 2a + 7 Li 8.82 

n + 16 — > n + 4a 14.44 

— > a + 13 C 2.22 

— > 2a + 9 Be 12.86 

— > 7 Li + 10 B 23.63 

n + 20 Ne — > a+ 17 0.59 

— > 2a + 13 C 6.98 

— > 3a + 9 Be 17.60 

— > 7 Li +a + 10 B 28.36 

— > 7 Li + 14 N 16.74 



